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ABSTRACT 



We present the first sky maps from tlie BEAST (Background Emission 
Anisotropy Scanning Telescope) experiment. BEAST consists of a 2.2 meter 
off axis Gregorian telescope fed by a cryogenic millimeter wavelength focal plane 
currently consisting of 6 Q band (40 GHz) and 2 Ka band (30 GHz) scalar feed 
horns feeding cryogenic HEMT amplifiers. Data were collected from two balloon- 
borne flights in 2000, followed by a lengthy ground observing campaign from the 
3.8 Km altitude University of Cahfornia White Mountain Research Station. This 
paper reports the initial results from the ground based observations. The instru- 
ment produced an annular map covering the sky from 33° < S < 42°. The maps 
cover an area of 2470 square degrees with an effective resolution of 23 arcminutes 
FWHM at 40 GHz and 30 arcminutes at 30 GHz. The map RMS (smoothed 
to 30 arcminutes and excluding galactic foregrounds) is 54 ± at 40 GHz. 
Comparison with the instrument noise gives a cosmic signal RMS contribution of 
28±3//K. An estimate of the actual CMB sky signal requires taking into account 
the £ space filter function of our experiment and analysis techniques, carried out 
in a companion paper (O'Dwyer et al. 2003). In addition to the robust detection 
of CMB anisotropies, we find a strong correlation between small portions of our 
maps and features in recent Ha maps (Finkbeiner 2003). In this work we de- 
scribe the data set and analysis techniques leading to the maps, including data 
selection, filtering, pointing reconstruction, mapmaking algorithms and system- 
atic effects. A detailed description of the experiment appears in Childers et al. 
(2003). 

Subject headings: cosmology: experiments 
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1. Introduction 

The Cosmic Microwave Background (CMB), interpreted as relic radiation from an 
early hot dense phase of the universe, is generally acknowledged as one of our primary 
probes of conditions in the early universe(Hu and Dodelson 2002 and references therein). 
Mapping the temperature anisotropics of the CMB has been the subject of intensive effort 
over the past two decades, with the results of many studies published within the last 3 
years: BOOMERANG (Netterfield et al. 2001), MAXIMA(Lee et al. 2001), CBI (Mason ct 
al 2002), DASI(Halverson et al. 2002), Archeops (Benoit et al. 2002), and ACBAR(Kuo et 
al. 2002) among others (see BersaneUi et al 2002 for a recent review). 

Temperature anisotropy measurements have been used to calibrate the amplitude of 
the initial density perturbations and constrain the primordial spectral index (Gorski et 
al. 1996), measure the overall density of the universe(De Bernardis et al. 2001, etc.), and 
recently to make detailed comparisons between competing scenarios for structure formation, 
leaving a rather narrow range of consistent adiabatic inflationary models. In combination 
with recent constraints from supernova SNIa and Large Scale Structure (LSS) data, CMB 
temperature anisotropy data strongly suggest models with a large fraction of the total 
energy density in a cosmological constant or other 'dark sector'. 

CMB data are normally reduced to an estimate of the angular power spectrum of 
fluctuations for comparison to theoretical predictions. The salient features of the currently 
allowed adiabatic inflationary models are a series of 'acoustic' peaks and troughs in the 
power spectrum from £ (spherical harmonic multipole) of about 200 and up. (Bond & 
Efstathiou 1987) The positions and amplitudes of these peaks pin down the details of the 
cosmological models. Each of the aforementioned experiments was designed to target a 
particular region of '^-space' or portion of the angular power spectrum. Sensitivity to high 
I modes is limited by the angular resolution of the telescope. Resolution in £-space is 
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limited by the size of the region studied. Amphtude uncertainties at all I are limited by 
the detector noise and integration time. At low £ values in particular, sample variance can 
increase the amplitude uncertainty of a measurement, an effect which can only be reduced 
by increasing sky coverage. BEAST, in particular the observing campaign described in this 
paper, has been optimized to study the first of these peaks in detail. 



1.1. Instrument Overview 

BEAST is a mixed Q-band (38-45 GHz) and Ka-band (25-35 GHz) focal plane that 
views the sky through a 2.2 m off-axis telescope followed by a modulating flat mirror. 
Primarily funded by NASA, BEAST is largest aperture CMB telescope flown to date. 
The instrument is described in detail in Childers et al. (2003), while the optical design 
is described in Figueiredo et al. (2003). BEAST is most similar to one of our previous 
balloon borne telescopes, HACME (Staren et al. 2000; Tegmark et al. 2000) but with a 
much larger telescope and a focal plane array. After two balloon flights in 2000, BEAST 
was reconflgured to take advantage of the UC White Mountain Research Station (WMRS) 
site on Barcroft mountain, CA (altitude 3800 m ). The data described here were taken with 
6 elements of the array operational, 2 Ka-band and 4 Q-band detectors. 2 of the Q-band 
detectors are installed on an Ortho-Mode Transducer (OMT) using a single horn split into 
vertical and horizontal polarizations. The other horns have single polarization receivers, 
oriented at various angles. The telescope was installed in an unused building, modifled so 
the roof can be rolled off. BEAST is ahgned at 74 degrees East of North, and can view 
nominal central elevation angles from 60 to greater than 90 degrees. The instrument was 
fully installed and operational at Barcroft in July, 2001, and took data nearly continuously 
until December 2001. Two more weeks of data were obtained in February 2002 and a longer 
campaign took place from August through October of 2002. 
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1.2. Observing strategy 

BEAST views the sky through a large (2.5 meter diameter) spinning flat mirror, tilted 
2.2 degrees from its rotation axis. As the flat spins this tilt moves the field of view of the 
telescope around a nearly elliptical path on the sky of diameter about 9 degrees. The 5 
operational feeds thus follow the path shown in Figure 1 for the normal telescope elevation 
angle of 90.4 degrees. Two main strategies were employed in taking data. The first involved 
scanning the nominal (central) elevation of the telescope from 80 to 90 degrees and back 
continuously, the second method was to keep the elevation angle fixed near 90 degrees and 
allow Earth rotation to provide the mapping scan. All of the data reported here are from 
the 'fixed elevation' data set. This choice was made primarily for the simpler and more 
robust pointing solution obtained in fixed elevation mode, and the fact that there were no 
significant improvements in systematic tests gained from the elevation scans. The resulting 
map is an annulus around the North Celestial Pole (NCP) of width ~ 10 degrees, around 
declination 37 degrees. Figure 2 shows the entire Q-band map with the beam trajectories 
overlaid. 

The output of each radiometer is amplified, ac coupled with a highpass filter time 
constant of 15 seconds, and then integrated in an 'ideal' or boxcar integrator with nearly 
100% duty cycle. The integration time is set for 800/xs and digitized to 16 bits. 

The angular position of the flat mirror (rotating at 2 Hz) is measured simultaneously 
so the data can be synchronized with the instantaneous position of the relevant telescope 
beam. The first step of analysis involves rebinning the raw 1250 Hz data stream into 
'optical sectors', corresponding to 250 angular positions of the flat mirror. Each 'optical 
sector' therefore, identifies an angular position of the fiat, which is related by the telescope 
geometry to a precise azimuth-elevation pointing for each of the 5 active detector elements. 
These pointing data are converted to celestial coordinates and then ultimately to HEALPix 
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(Hierarchical Equal Area isoLatitude Pixelisation, Gorski et al. (1998)) pixel numbers. The 
resulting optical sector data stream has an effective sampling rate of approximately 450 Hz. 
The raw data amount to approximately 3.7 Gigabytes per day. 

1.3. Calibration 

The receivers are calibrated hourly by passing an ambient temperature microwave 
absorber between the feeds and the secondary mirror. The output of the acquisition 
system is then corrected for the effects of the ac coupling in addition to a small (order 
10%) correction for compression of the detector outputs while viewing the ambient target. 
The difference between the load and the sky temperatures is then used to compute the 
calibration coefficients, measured in K/V. This hourly calibration constant is applied 
to the data during pipeline analysis. Calibration and effective beamsize were confirmed 
independently using Cygnus A as the celestial source. Overall calibration uncertainty is 
±10%. Calibration and beamsize issues are covered more completely in Childers et al. 
(2003). 

2. Data Set 

2.1. Raw Data Statistics 

Figure 3 shows a sample section of time ordered data for a typical channel. The top 
strip of small points is the raw data stream, including electrical and any radiometric offset. 
The large nearly sinusoidal component is due to the varying atmospheric contribution as 
the fiat sends the beam through different elevations and therefore different optical depths. 
This variation is approximately 60 mK peak-to-peak for this Q-band (38-45 GHz) channel. 
We obtain a high signal to noise estimate of the nearly constant atmospheric contribution 
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for each revolution by binning all the data for each hour by optical sector to produce a 
'template' signal for that hour. We then subtract this template from every revolution in 
that hour to remove the first order atmospheric contribution. Note that systematic errors 
fixed in the telescope frame or phase locked to the flat revolution would be removed as well. 
Examples of such effects are flxed sidelobe contamination and electrical noise synchronous 
with the flat motion. Since our strategy keeps the telescope flxed, in principle any sidelobe 
illumination of the observatory or ground would be removed. There is no direct evidence for 
any contribution except atmosphere. The solid line plotted through the raw data points is 
the template for this hour, while the central strip of data points clustered near zero are the 
residuals from the template subtraction. The lower strip of data points show the result of 
highpass flltering the template-removed data set at 10 Hz corresponding to about 6 degrees 
on the sky (effectively giving us a low £-space cutoff) . These data have been offset from 
zero for clarity. The combination of 1 hour template removal and steep 10 Hz highpass fllter 
is a reasonable compromise between between extending the low £ coverage and minmizing 
the increasing effects of 1/f instrument and atmospheric noise contributions. Other fllters 
have also been tried and will be reported on in the future. Without the template removal, 
harmonics of the spin rate remain in the data and can show up as striping in the maps. 
Hourly sections of data are our natural units since calibrations happen on the hour, breaking 
the continuity of the data stream. There appears to be no beneflt to using longer sections 
for calculating the template: we have not yet investigated in detail the advantages of using 
shorter time periods for template subtraction. Both hourly template removal and 10 Hz 
highpass fllter were used in the analysis described below. 

Estimates for the noise amplitude spectral density for this hour are given in Figure 
4, cahbrated in jjK^Js. The large signal near 2 Hz and its higher frequency harmonics 
are due to the atmospheric contribution discussed above. Notice that due to fact that 
the beam trajectory is not circular, and passes beyond zenith, the resulting atmospheric 
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Tablc 1: Summary statistics for a typical hour of sky data 



Channel 


Frequency 


White Noise 


1/f Knee 


Template p-p 




(GHz) 


(//Kv^) 


/tt \ 

(Hz) 


(mK) 


2 


38-45 


697 


38 


84 


3 


38-45 


1324 


40 


149 


4 


38-45 


770 


60 


102 


6 


25-35 


1106 


105 


76 


7 


25-35 


966 


110 


107 


8 


38-45 


642 


58 


87 



signal is distorted from sinusoidal, producing significant harmonics of the spin frequency. 
Also shown on the plot is the noise density after template removal and 10 Hz highpass 
filtering. The plot demonstrates the BEAST detector noise has significant contributions 
from correlated, long memory (or "1/f") noise. The knee frequency of "1/f" noise is defined 
as the frequency where the total noise power spectral density (1/f + white) exceeds by a 
factor of two its white noise level. The latter determines our sensitivity per unit time for 
angular scales near the beam size. Table 1 summarizes the receiver and atmospheric noise 
and signal performance for a typical accepted hour of data from this campaign. 



2.2. Pointing Reconstruction 

The BEAST scan strategy produces a large area map (2470 square degrees) with high 
resolution, and we require a precise pointing solution to associate every optical sector 
for every channel with a position in celestial coordinates. We know the geometry of the 
telescope, but not well enough to calculate the beam positions to sufficient accuracy for 
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map reconstruction. Deviation from ideality of the optics gives uncertanties on the "a 
priori" knowledge of the beam positions. BEAST was placed in an existing building that 
did not allow us to move the telescope appreciably in azimuth, and we are unable to target 
the moon, sun or any planets for pointing verification. We bootstrapped the pointing model 
by starting with the known geometry, then flying a small airplane equipped with GPS over 
the system in a grid pattern. The spikes in the BEAST data due to the airplane thermal 
emission were then correlated with the measured relative position of the airplane to make a 
crude improvement to the pointing model. Cygnus A as well as several other well known 
radio sources are included in our map, and we have been able to use these sources to 
further refine the pointing model, beamsize and calibration. Figure 5 is a portion of the 
final map (all Q-band channels averaged) showing Cygnus A and the Cygnus Loop region. 
The reconstruction returns an effective resolution of 23 arcminutes FWHM in Q-band. This 
is in contrast with the previously measured telescope resolution of 19 arcminutes. The 
smearing from the design resolution to the measured effective FWHM of 23 arcminutes is 
due to a combination of unmeasured pointing errors (telescope fiex, long term telescope 
sag), residual errors in the pointing reconstruction algorithm, smearing due to the finite 
HEALPix resolution of 6.9 arcminutes, and smearing from the inital flat rotation sectors 
(about 6.7 arcminutes). The pointing model is converted to a lookup table of azimuth and 
elevation for each horn, sector and telescope elevation for use in the pipeline processing. 

2.3. Data Reduction and Analysis 

The first data reduction step was to reformat the data into optical sectors, and 
incorporate relevant pointing information into a single file (hereafter called 'level FITS 
file'), which conforms to the FITSIO standard. A single file was produced for each hour, 
including all channels, as well as telescope elevation, and other pointing information. In 
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addition to rebinning to optical sectors we measured some statistics on these hourly sets to 
allow simple data culling later. In particular we measured the RMS of each channel for the 
hour as well as the RMS after highpass filtering the data at 10 Hz. Transients (generally 
due to thermal emission from airplanes crossing the instruments field of view) were removed 
and replaced with adjacent data. These points were flagged as 'unusable' for later analysis 
stages. Level processing was:l) Import 1 hour Receiver and Pointing data, 2) Calculate 
statistics (RMS and highpass filtered RMS), 3) Rebin raw (not highpass filtered) data to 
Optical sectors, 4) Remove transients, flag bad or replaced points 5) Store data to FITS, 
save summary statistics in header. 

The next step was to cull the data for bad weather and out of specification performance 
(ie when systems were off or under test or during calibrations). This was done by looking 
at the statistics of the highpass filtered data RMS in the level headers. Generally all the 
channels showed the same sensitivity to bad weather, and there was a clear knee in the 
distribution of RMS's for making the cut. 

We then moved to 'level 1' data files, which integrated the pointing data into map 
pixel indices. We chose to use the HEALPix (Gorski et al. 1998) set of standards and 
routines, partly for the inherent advantages of the pixelization scheme and partly for the 
extensive infrastructure of tools available from the HEALPix team for manipulating and 
visualizing the resulting maps. This standardization and set of tools has been invaluable 
in performing all of our subsequent analysis Every sample for each channel was then 
associated with an azimuth-elevation pointing via the lookup table mentioned above, then 
to a celestial coordinate pointing using the known observation time and location, thence to 
a HEALPix index using the HEALPix IDL library routines. We chose to use Nside=512, 
with a resolution corresponding to 6.9 arcminute pixels as a reasonable oversampling 

-^^The HEALPix homepage is located at http://www.eso.org/science/healpix/. 
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of our beamsize. Level 1 data files then included time, 6 columns of data, 1 column of 
'goodness-of-data' flag, and 6 columns of HEALPix indices. 

To produce a map, we 1) Read in Level 1 file, 2) Calculate and subtract hourly 
atmospheric template from all channels, 3) highpass filter all channels at 10 Hz, 4) Bin into 
daily map (ignoring flagged points, saving statistics), 5) Iterate for all hours to complete 
daily map for all channels, 6) Save daily map to disk, start next day This generates a map 
with statistics for each day for each channel. We then co-added (averaged by pixel) the 
daily channel maps into total data set channel maps, weighting by the average noise per 
pixel per day. The Q-band channels were then co-added with appropriate noise weighting to 
create the total Q-band data maps, as were the Ka channels. In addition to these averaged 
maps, we made 'difference' maps, where we binned 1/2 of the days into one map, the other 
1/2 of the days into another. Then we generated the difference of these two maps on a pixel 
by pixel basis, dividing by 2 to maintain the noise statistics of the average maps. These 
difference maps give us a measure of the noise in the maps vs the signal -|- noise in the 
average maps. Some of these difference maps are displayed below. 



2.4. Maps and Statistics 

Figures 6 and 7 display two representative sections of the BEAST Q-band average 
map. Each Figure contains 3 panels. The top panel is the full Nside=512 resolution map in 
gnomonic projection around 5 — 3S degrees and a— and 160 degrees. The second panel 
is this same map smoothed with a 0.5 degree FWHM gaussian. The third panel is the 0.5 
degree FWHM smoothed difference map as a visual guide to the effective noise level in the 
smoothed average map. 



2.5. Noise properties 



Figure 8 shows the distribution of pixel values for the Q-band map, compared with a 
gaussian distribution. The distributions with and without the galaxy and point source cuts 
are shown. One measure of the noise and signal in the map is to calculate the standard 
deviation of the map (after removing foreground contaminated regions) and compare 
with the standard deviation in the difference map, which should have the same noise 
characteristics, but no CMB signal. For Q-band we find 254 ± 25/xK at 6.9' pixels for the 
sum map and 252 ± 25/xK for the difference map. Smoothing to 30 arcminutes FWHM 
we get 54 ± 5//K for the sum and 47 ± 5//K for the difference, consistent with a signal 
RMS contribution of 26/x K (Rayleigh- Jeans units). Converting to blackbody spectrum and 
correcting for 3% atmospheric attenuation we obtain 28 ± 3/xK (quoted error dominted 
by calibration uncertainty). Note that this includes the map smoothing, and the complex 
BEAST filter produced by our AC couphng, beam shape, highpass filter, and template 
removal, so is not our best estimate of sky signal, but only an indicator of a statistically 
significant sky signal. 

2.6. Foregrounds 

Any sensitve CMB map of this size will cross regions of significant foreground 
contamination as is seen in our maps with the galactic crossings. Other than a few regions 
the map appears to be remarkably clean. By comparison with foreground maps, we have 
been able to determine a very conservative cut in galactic latitude which removes diffuse 
foregrounds, while for point sources we have used a strict point source finding algorithm 
operating on the data maps. Point sources were removed by excising a 1 degree diameter 
circle of pixels centered on each source. The most interesting feature of our foreground 
contamination is a large structure near a— 60 degrees, 5=36.5 degrees. The feature is 0.75 
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mK in the smoothed Q-band map and 1.4 mK in Ka-band, implying a power law spectral 
index of -1.9, consistent with free-free emission. We see a large correlation with a feature 
in the combined IRAS/DIRBE map (Finkbeiner et al. ), but no correlation at all with 
the Haslam map at 408 MHz. Initially we thought this might be evidence for spinning 
dust grains, but a comparison with integrated Ha emission showed a very large correlation, 
consistent with free-free emission. This comparison was made possible by the availability 
of a map of integrated Ha emission produced by D. Finkbeiner (Finkbeiner, 2003) using 
data from WHAM (Reynolds et al. 2002), VTSS (Dennison et al. 1998), and SHASSA 
(Gaustad et al. 2001). The expected free-free emission is estimated by scaling from Bennett 
et al. (1992) which predicts 3/iK /Rayleigh at 40 GHz. For our structure we measure 
5.7/iK/Rayleigh. Since the correlations with the Ha map were so much more evident than 
with Haslam and IRAS that we used the Ha data to estimate our required galactic latitude 
cut: by cutting out data with | b \ less than 25° we expect less than 10/xK contamination 
from free-free emission. This is a very conservative cut and we expect to refine this further 
to increase the accepted sky fraction. Our resulting foreground removal template is shown 
in figure 9, shaded regions are used in the power spectrum analysis. 

3. Discussion 

We have produced extended (~ 2500 square degree) sky maps with good angular 
resolution and control of systematics. Visual inspection of the smoothed maps shows 
unmistakable structure. Quantitative results for the power spectrum of CMB fluctuations 
which is consistent with our data are presented in a companion paper (O'Dwyer et al. 2003). 
In addition to cosmic structures we have interesting correlations with known astrophysical 
foregrounds. Further analysis will yield important information about the amplitude, spectra 
and possibly polarization of these foregrounds. 
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As an aid in interpreting our maps, and as a guide to our ^-space filter, we have made 
a single realization of a 'white' power spectrum sky (equal power on all scales), and then 
simulated the BEAST pipeline over that map. We created a full sky map with normally 
distributed noise in each pixel. This map was then smoothed with a gaussian of FWHM=23 
arcminutes to simulate our effective beam. Next we simulated the time ordered data 
expected from this map without additional noise by sampling with the BEAST pointing 
timestream. We ran the resulting timeline through the BEAST analysis pipeline, including 
template removal, filtering and binning. Finally we cut using the same foreground template 
used in the actual data map. 'PseudoCl' power spectra (that is power spectra calculated 
assuming full sky coverage but actually covering only a particular region) were calculated 
for the full unsmoothed white input map, the smoothed white map, the smoothed and 
foreground cut white map, and the full BEAST simulated smoothed and foreground cut 
map. The results are shown in figure 10. The two all sky power spectra have been scaled 
by the sky fraction (2.85% of full sky) of the foreground cut template in order to put them 
onscale, and to demonstrate that the primary effect of limited sky coverage on our £-space 
filter is just amplitude. All of the smoothed maps display the high £ cutoff of our beam. 
Below £ of about 150, the input map power spectrum is virtually the same as the smoothed 
and cut sky power spectra, while the BEAST simulated map power spectrum displays 
the low i cutoff of the combined template removal and 10 Hz highpass filter. The loss of 
sensitivity from £ ?a 50 to £ ~ 150 is dominated by the highpass filter cutoff. Note that this 
analysis neglects correlations among different £s. Figure 11 gives an approximation to the i 
space window function based on these calculation. 

Power spectra from the data and cosmological parameter estimates will be presented 
in a companion paper (O'Dwyer et al. 2003). 
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4. Future Observations 

We have shown the utihty of making large area sky maps with a simple scanning 
strategy and expect to add significantly more data in the near future. The data presented 
so far are the result of 27 days of "good data". Our experience so far is we can achieve 
"good data" more than 50% of the days at our observatory. The site has proven to be 
extremely good in our wavelength range and the atmosphere stable. A paper discussing 
the site is in preparation. One year of observations with the existing system only improved 
by new HEMT's should yield 7 /jK errors per 30 arcminute pixel on 10% of the sky. The 
BEAST optics are precise enough to observe to nearly 1 THz so a variety of sky surveys are 
possible. At 90 GHz the beam size would be about 8.5 arc minutes FWHM while at 150 
GHz it is 5 arcminutes. 
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Fig. 1. — Trajectories of the BEAST main beams. The figure is a polar plot about the local 
zenith. We subdivide the path of each beam into 250 optical 'sectors' parameterized by the 
phase angle of the flat mirror which spins at about 2 Hz. 
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Fig. 2. — The BEAST combined Q-band map. The figure is a gnomonic projection about 
the north celestial pole. The grid lines are 5 degrees in declination and 5 degrees in right 
ascension. The trajectories of the 5 pixels during 1 revolution of the fiat mirror are overlaid 
on the figure near a=0° to illustrate the observing strategy. The fiat spins at about 2 Hz: 
the map includes data from over 4 million flat revolutions. Two galactic plane crossings are 
clearly visible on the plot, the brigher around a— 20°. Also visible is the large structure at 
a fa 60° found to be due to free free emission (see the text) . 
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Fig. 3. — Sample of time ordered data (a; 8s). The top strip of small points is the raw 
data stream. The solid line through the points is the flat synchronous atmospheric template 
obtained by averaging over this hour (see the text). The central strip of data points clustered 
near zero are the residuals from the template subtraction. The lower strip of data points 
show the result of highpass filtering the template-removed data set at 10 Hz, offset from zero 
for clarity. 
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Fig. 4. — Noise amplitude density spectrum from the sample hour of data. The white noise 
limit at high frequency, the large 1/f component, the atmospheric signal at the spin rate 
and its harmonics, and excess noise at 120 Hz are all evident. The noise spectrum for the 
template subtracted data (dotted line), and the spectrum after 10 Hz high are plotted as 
well(dashed line). 
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Fig. 5. — Gnomonic projection of the Q-band map near Cygnus A and the Cygnus Loop 
region. The grid spacing is 5 degrees in RA and Dec. Cygnus A has been used as a pointing 
cahbration as well as a cross check of the effective beamsize and calibration constant. All of 
the sources seen here have been identified with known objects. 




Fig. 6. — Gnomonic projection of the Q-band map centered at a= 0° and 6= +38°. The 
top map is the full resolution (6.9 arcminute pixels) map. The center map is smoothed with 
a 30 arcminute FWHM gaussian. The bottom map is the first half/second half difference 
map, also smoothed to 30 arcminutes. The difference map should have only noise while the 
average map (center) should show noise+ cosmic structure. Declination increases up, right 
ascension increases to the left, the grid is 5 degrees spacing in right ascension and 10 degrees 




Fig. 7. — Gnomonic projection of the Q-band map centered at a= 160° and 6= +38°. 
This section of the map has an obvious point source at the upper right (removed for power 
spectrum analysis) 
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Fig. 8. — Normalized histogram of the full resolution Q-band map pixels. Symbols are the 
full map, the solid black line is the map with galaxy and point sources removed (| & |< 25°), 
the dashed line is a gaussian with 250/xK a. 
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Fig. 9. — The map sections available for power spectrum analysis after foreground cuts 
(I ^ l< 25°) 
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Fig. 10. — An approximate ^ -space description of the effective filtering produced by our 
pipeline processing and galaxy cuts. Starting at the top: Power spectrum for the full 'white' 
simulated sky, scaled by the sky fraction; Spectrum for 23 arcminute FWHM smoothed 
simulation, scaled by sky fraction; 23 arcminute smoothed simulation cut by galaxy removal 
and BEAST sky coverage template; Smoothed sky simulation measured through BEAST 
pipeline analysis (no noise added). Note that all spectra were caluclated using the HEALPix 
'ANAFAST' routine. Also all except the full sky spectrum are 'pseudo-Cl' spectra, calculated 
from partial sky coverage without any corrections. See the text for discussion. 
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Fig. 11. — The effective I space window function of the experiment and analysis. Plotted 
is the relative sensitivity of the BEAST experiment including beam smoothing, highpass 
filtering and template removal, as compared with a perfect experiment measuring the same 
region. The comparison is done using the same methods as the previous figure. This cal- 
culation does not include the effects of the noise, only the instrument and analysis filtering. 



